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The  Effect  of  the  Convergence  of  the  Meridians  on  the 
Westward  Intensification  of  the  Oceanic  Circulation 

Abstract: 

^Models  for  a wind  driven  ocean  were  studied  including  the  effects 
of  sphericity  and  rotation  of  the  earth.  Like  Stommel's  1948  model,  they 
exhibited  the  westward  intensification  of  the  circulation,  but  significant 
difference  in  pattern  appeared. 

Introduction: 

Stommel  (1948)  established  that  the  variation  of  the  Coriolis  para- 
meter with  latitude  is  responsible  for  the  westward  intensification  in  an 
ocean  of  homogeneous  water.  In  that  model,  he  assumed  that  the  varia- 
tion is  uniform  with  latitude.  A study  was  initiated  to  investigate  the  ef- 
fects of  dropping  that  assumption.  Further,  the  inclusion  of  the  Coriolis 
parameter  variation,  makes  rectangular  ocean  models  bounded  by  meri- 
dians and  latitude  circles  on  a rotating  earth  inconsistent,  as  it  can  only 
change  where  meridians  converge.  Two  further  models  were  therefore 
considered,  where  the  meridians  converge.  One  modelled  a conical  earth, 
consistent  with  a unform  change  in  the  Coriolis  parameter,  and  the  other 
modelled  a spherical  earth. 

Equations : 

The  treatment  was  kept  as  near  to  Stommel's  1948  treatment  as 

possible,  and  most  of  the  original  symbols  are  thus  retained.  Assume 

an  ocean  bounded  by  two  meridians  and  two  latitude  circules  with  corners 

at  (0,0),  (A.0),  (X,b),  (0,b).  Assume  further  that  a zonal  wind  stress 

steadily  blows  over  the  ocean  and  is  expressed  by  T = - F cos-^  . 

x b 

The  equations  of  motion  thus  are: 
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f (D  + h)  v p - F cos  (TTy  /b)  - K u (D  + h)-  g p(D  + h)  ^ = 0 


(1) 


-f(D+h)up  - Kv(D+  h)-gp(D+  h)  |^=  0 (2) 

The  equation  of  continuity  is: 

(a)  for  the  case  where  the  convergence  of  the  meridians  is  ignored 

J^((D  + h)  u)  + f^((D  + h)  v)  = 0 (3a) 

(b)  for  the  conical  earth 

f^0D  + h)  0 + ^fy  0D  + h)  v)  = ~y^yQr  (D  + h0  (3b) 

(c)  for  the  case  of  a spherical  earth 

f^(u  (D  + h))  + (v  (D  + h))  = -*522(v  (D  + h))  (3c) 


Equations  3a,  3b,  3c  allow  the  introduction  of  a stream  function  ^ such 
that 


u (D  + h)  = - , v (D  + h)  = + 

(4) 

j 

Utilizing  this  relation  a vorticity  equation  can  be  developed 

K72t>pp|4  ^ 

(5a) 

(5b) 

vr  n 2^  J.  20  pcos  (P  Ktan<pj£&_  0Zx  + tan^j 

K R 8*  ■ R »y‘  8y  R T* 

(5c) 
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the  length  of  the  ocean  along  a meridian 
the  Coriolis  parameter 
the  acceleration  of  gravity 

the  displacement  of  the  sea  surface  from  its  undisturbed  level 

the  velocity  in  the  zonal  direction 

the  velocity  in  the  meridional  direction 

the  zonal  coordinate,  increasing  eastward 

the  meridional  coordinate,  increasing  northward 

the  coordinate  of  the  pole  in  the  conical  earth 

the  depth  of  the  ocean 

the  maximum  wind  stress 

the  coefficient  of  friction 

the  radius  of  the  earth 

the  variation  of  Coriolis  parameter  with  latitude 

+ J&L 

9x2  8y2 

the  maximum  width  of  the  ocean  considered 

the  density  of  the  water 

the  wind  stress  s F cos  ^ 

b 

the  latitude 

the  transport  stream  function 


3 


Equations  (5a),  (5b),  (5c)  were  solved  numerically  and  the  solu- 
tions appear  in  Figure  1,  2,  3,  4.  Figure  1 reproduces  Stommel's  pic- 
ture, while  2 is  for  the  same  model  with  a variable  3.  Figures  2 and  4 
are  for  a conical  and  for  a spherical  earth  respectively.  Dimensions  of 
the  oceans  and  relevant  parameters  appear  in  Table  I.  Latitudes  are 
given  on  all  of  the  figures,  though,  they  are  relevant  only  in  Figures  2 
and  4. 

Table  I:  Dimensions  of  oceans  and  other  relevant  parameters. 

o 

b = length  of  ocean  along  a meridian  = 0.5  X 10  cm 

9 

X = maximum  width  of  ocean  =10  cm 

9 

yQ  - the  coordinate  of  the  pole  in  the  conical  earth  = 10  cm 

D = the  depth  of  the  ocean  = 200  m 

F = maximum  wind  stress  = 1.0  dyne  cm 

K = coefficient  of  friction  = 10"^  dynes  cm  4 sec 

9 

R = radius  of  earth  = 0.6371  X 10  cm 

(3  = (when  used  as  a constant)  = lO’*"*  cm4  sec"* 

-2 

p = density  of  water  = 1 gm  cm 
Discussion  of  the  Results; 

The  circulation  remains  intensified  in  the  west,  but  the  pattern 
shows  interesting  changes.  When  the  convergence  of  the  meridians  were 
ignored,  the  center  of  the  gyre  was  displaced  to  the  north.  The  transport 
decreased  as  compared  with  the  case  of  a constant  3.  It  is  to  be  noted, 
however,  that  the  mean  (when  3 varied)  was  greater  than  the  constant 
value  that  was  used  corroborating  the  results  of  a study  made  by  Garner 
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(1962)  about  the  relation  between  the  value  of  (3  and  the  transport. 
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When  the  convergence  of  the  meridians  were  taken  into  account, 
however,  the  center  of  the  transport  gyre  was  displaced  to  the  south. 
The  value  of  the  transport  was  again  decreased,  reflecting,  in  addition 
to  the  (3  effect,  the  effect  of  the  decrease  in  area  affected  by  the  wind. 
A weak  circulation  with  the  opposition  sense  of  rotation  also  appeared 
at  the  northwestern  corner  of  the  model. 
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Second  modified  model,  conical  earth,  constant  (3.  Stream  functions 


Figure  4.  Third  modified  model,  spherical  earth,  variable  p.  Stream  functions 
in  c. g. s.  X 10+9. 
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Approximations  of  the  Nonlinear  Terms  in  the  Vertically 
Integrated  Equation  of  Motion 

Abstract: 

Integrated  nonlinear  terms  in  the  hydrodynamical  equations  of 
motion  were  compared  with  the  nonlinear  terms  formed  from  the  inte- 
grated velocity.  This  showed  a difference  in  form  to  indicate  an  essen- 
tial difference  between  the  two  expressions.  Numerical  values  of  the  two 
expressions  were  compared  for  a particular  model  and  the  results  sub- 
stantiated the  theoretical  treatment. 

Introduction: 


The  nonlinear  terms  in  the  hydrodynamical  equations  of  motion 
pose  a challenge  when  the  equations  are  considered  for  solution.  For  a 
variety  of  reasons,  the  equations  become  more  amenable  to  treatment 
when  they  are  integrated  vertically.  This  poses  an  additional  difficulty 
as  there  is  no  simple  expression  relating  the  vertically  integrated  non- 
linear terms  with  those  which  can  be  formed  from  the  vertically  integrated 
velocities.  It  has  been  occasionally  assumed,  however,  (Bryan  1963, 
Carrier  and  Robinson  1962)  that  a simple  proportionality  exists.  This  is 
investigated  here. 

The  velocity  is  expanded  in  the  vertical  in  a Fourier  series  in 
such  a way  that  the  velocity  at  the  bottom  automatically  vanishes  without 
imposing  any  restrictions  on  the  velocity  profile.  Expressions  are  then 
evaluated  for  the  integrated  nonlinear  terms  and  the  nonlinear  terms 
formed  from  the  integrated  velocities.  These  terms  are  then  compared. 
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Procedure: 


Assume  that 


a=I  “i 


( ii  - 1 )V  z 


Then 


and  u 

dx 
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It  is  obvious  that  (1)  and  (2)  do  not  have  a simple  proportionality 
factor.  To  get  a further  feeling  for  these  terms,  an  oceanic  model  was 
used,  for  which  an  analytic  solution  could  be  produced  (Hassan,  1963). 

The  model  has  many  deficiencies,  the  most  relevant  here  is  that  it  does 
not  contain  the  nonlinear  terms.  The  terms  computed  were  the  following: 

j*  fudz,  [ fv  dz,  |u|j  dz,  j~u  dz,  Jv  §£  dz,  §£  dz, 

o o o o o o 

h h h h h h 

h £ddI  fe  h iz-  h d*’ 

o o o o o o 

h h 

h dz  % Jvdd- 

o o 

The  first  two  terms  were  computed  to  indicate  whether  the  model 
can  still  be  used  to  evaluate  the  nonlinear  terms,  even  though  the  non- 
linear terms  were  not  considered  in  the  solution.  Not  surprisingly,  ex- 
cept for  the  extreme  western  boundary,  the  nonlinear  terms  were  negli- 
gible compared  to  the  terms  retained,  and  only  at  the  western  boundary, 
were  they  of  comparable  magnitude  (Figure  1 and  2).  The  two  expressions 
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= I 

for  the  nonlinear  terms  showed  very  little  relationship  to  one  another, 
as  can  be  seen  from  Figures  3 to  10. 

Conclusion: 

Vertically  integrated  nonlinear  terms  cannot  be  approximated  by 
expressions  using  the  vertically  integrated  velocities. 
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